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Abstract 

The geodesic deviation equation ('GDE') provides an elegant tool to investigate the timelike, null 
and spacelike structure of spacetime geometries. Here we employ the GDE to review these structures 
within the Friedmann-Lemaitre-Robertson-Walker ('FLRW') models, where we assume the sources 
to be given by a non-interacting mixture of incoherent matter and radiation, and we also take a 
non-zero cosmological constant into account. For each causal case we present examples of solutions 
to the GDE and we discuss the interpretation of the related first integrals. The de Sitter spacetime 
geometry is treated separately. 
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1 Introduction 

It has been known for a long while that the geodesic deviation equation ('GDE'), first obtained by J L 
Synge ||2^, 24], provides a very elegant way of understanding features of curved spaces, and, as pointed out 
by Pirani [14, 15 , gives an invariant way of characterising the nature of gravitational forces in spacetime. 
As such, it is a useful tool to use in examining specific exact solutions of the Einstein field equations 
('EFE'). Indeed, it may be claimed that the GDE is one of the most important equations in relativity, 
as this is how one measures spacetime curvature^. This latter aspect has been discussed in some depth 
by Szekeres |^ . 

The GDE determines the second rate of change of the deviation vectors for a congruence of geodesies 
of arbitrary causal character, i.e., their relative acceleration. Consider the normalised tangent vector field 
V"" for such a congruence, parametrised by an affine parameter v. Then 

V^:=^^, VaV^:=e, 0=^-^^V'VtV\ (1) 
av ov 

where e = + 1, 0, — 1 if the geodesies are spacelike, null, or timelike, respectively, and we define covariant 
derivativion along the geodesies by 6T°"-b,J5v := V^V cT°'"b.. for any tensor T"- ;,... A deviation vector 
77" :— dx°'{w)/dw for the congruence, which can be thought of as linking pairs of neighbouring geodesies 
in the congruence, commutes^ with V"" , so 

^ = rfV^V- . (2) 
ov 

It follows that their scalar product is constant along the geodesies: 

5(n ¥"■) 

'- =0 ^ iVaV)^ const . (3) 

ov 

To simplify the relevant equations, we always choose them orthogonal: 

r]aV'' = 0. (4) 

The general GDE takes the form 

^ = -R"bcdV'v''V'' , (5) 

(see, e.g., Synge and Schild Schouten or Wald ). The general solution to this second-order 
differential equation along any geodesic 7 will have two arbitrary constants (corresponding to the different 
congruences of geodesies that might have 7 as a member). There is a first integral along any geodesic 
that relates the connecting vectors for two different congruences which have one central geodesic curve 
(with affine parameter v) in common. This is 

?7ia-^ ~ri2a-^ = const , (6) 
and is completely independent of the curvature of the spacetime manifold. 



The aim of this paper is to systematically use the GDE to explore the geometry of the standard 
Friedmann-Lemaitre-Robertson-Walker ('FLRW') models of relativistic cosmology (see, e.g., Refs. [^8[ 
|30| , H ), solving the GDE for timelike, null and spacelike geodesic congruences in these geometries; hence, 
obtaining the Raychaudhuri equation |Q determining the time evolution of these models Q , the Mattig 
observational relations ||ll| underlying the interpretation of cosmological data and determining the 
nature of their spatial 3-geometry ^. Also, we identify in each case the first integral for the GDE and 
comment on its meaning, in the null case leading to the usual reciprocity theorem , and in the timelike 
case obtaining generic solutions of the GDE via this integral. Thus, our purpose is to characterise the 
major geometrical and physical features of these spacetimes by use of the GDE, hence showing the utility 
of this equation in obtaining all the essential geometrical and dynamical results of standard cosmology 
in a unified way. It is a pleasure to dedicate this paper to Engelbert Schiicking, who has made a major 
contribution to obtaining clarity and elegance in understanding many features of relativistic cosmology. 

^And so is analogous to the Lorentz force law in electrodynamics; cf. Misner, Thorne and Wheeler, Ch.3, Box 3.1 ]l2| . 
^The Lie derivative of rj" along the integral curves of V" is zero; see, e.g., Schouten [El|. 
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1.1 The cosmological context 

In the cosmological situation we consider, we assume the sources of the gravitational field to be a non- 
interacting mixture of incoherent matter and radiation, to each of which the phenomenological fluid 
description applies ( see, e.g., Refs. |^ and Q ). For completeness we also include a cosmological constant 
A. 

Notation used is as follows: is the normalised timelike tangent vector field (ua m° = —1) to 
the fundamental matter fluid flow, which is geodesic: = u^\/i,u°' := ii"". The integral curves of m° 
are parameterised by the proper time t of comoving fundamental observers. We use standard FLRW 
comoving coordinates: 

ds^ = -dt^ + a^{t)f^^{xP)dx>'dx'' , f^^dx^'dx" = dr^ + f{r){d9^ +sin^ 9 d(f>^) , (7) 
= {dtT = ^"o , (8) 

where a{t) denotes the time dependent scale factor, and the function^ /(r) relates to the intrinsic curvature 
of the spacelike 3-surfaces {t = const} orthogonal to By spatial homogeneity and isotropy, the 
covariant derivative of reduces to 

VaUb = lehab, e := DaU'' = 3 - . (9) 
^ a 

Here, hat is the standard orthogonal projection tensor 

hab^9ab + UaUb hafj = ga/3 , (10) 

Q is the fluid rate of expansion, and the spatial derivative operator (projected orthogonal to on all 
indices) is denoted by Da ( cf . Ref . (9) ) . It is a well-known consequence of Eq. (^ that FLRW spacetime 
geometries have vanishing Weyl curvature (cf. Refs. ||^ and |^ ), 

Cabcd = ; (11) 

the fluid matter flow neither generates tidal gravitational fields nor causes propagation of gravitational 
waves. 

2 The Riemann curvature tensor 

In order to determine the exphcit form of the GDE (|), we need the Riemann curvature tensor Rabcd- 
Because of Eq. (|ll|), Rabcd can be expressed purely in terms of the Ricci curvature tensor Rab, its trace 
R, and the metric: 

Rabcd = \ ( Rac gbd — Rad Qbc + Rbd Qac ^ Rbc Oad ) — ^R{ Qac Obd — Qad 9bc ) • (12) 

The EFE algebraically determine Rab from the matter tensor Tab'^ 

Rab^Tab~\Tgab + Agab R=-T + 4A. (13) 

When the matter takes a 'perfect fluid' form: 

Tab ^ iti+p) Ua Ub+pgab T = - (/i - 3p) , (14) 

(/i is the total energy density and p the isotropic pressure) , the Ricci tensor expression is 

Rab = {tJ-+p)uaUb + ^{^i-p + 2A)gab ^ R ^ {fj, - 3p) + i A . (15) 
Thus, from Eq. (|2|), the curvature tensor takes the form 
Rabcd = I (m + A) ( gac gbd - gad gbc ) 

+ i (/i + p) ( gac Ub Ud - gad Ub Uc + gbd Ua Uc - gbc UaUd). (16) 



^Determined later by use of the 3-D spatial GDE. 

^Geometrised units, characterised by c = 1 = jr? ^ are used throughout. 
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Then, for any normalised vector field V^: Va = e, by a straightforward contraction one obtains from 
Eq. ( p^ ) the source term in the GDE: 

+ \{^i+p)[ {YbV^f gac - 2 (Hu') M(a K) +euau,]. (17) 

We will also want the GDE in the spacelike 3-surfaces {t = const} orthogonal to which are 3-spaces 
of maximal symmetry. In the FLRW case, the Gaufi embedding equation provides the relation 

^Rabcd = {Rabcd)± — | 6^ ( hac hbd — had hbc ) (18) 

for the 3-D Riemann curvature. From Eq. (p^, which made use of the EFE, one has 

iRabcd)± = i (/i + A) ( hac hbd - had hbc ) , (19) 

so that Eq. ( p8| ) becomes 

^Rabcd = K{t) ( hac hbd - had hbc ) , (20) 
where the spatial curvature scalar K(t) is given by 

K{t)■.= \^R=\{^^-\Q^ + k) . (21) 

This factor will determine the 3-D spatial GDE|^ source term: 

^Rabcd V' V^^K{ hac -VaVc) , (22) 

where VaV = 1 and Va = 0. 

3 The geodesies 

Before turning to address the GDE, we need to solve for the geodesic curves along which the GDE will be 
integrated. Now the fundamental 4-velocity u"- = 5% is a geodesic vector field. Any other geodesic can 
be transformed to have a purely radial spatial part by suitable choice of local coordinates (because the 
FLRW geometry is isotropic about every point). Hence, w.o.l.g., radial geodesies are considered, with 
the origin of the local coordinates r = at the starting point w = 0, so that in all cases we will have 
x'^ — 9 — const, = (f) — const => = = V'^ . 

It is convenient to decompose a general geodesic tangent vector field V"" into parts parallel and 
orthogonal to u"-; 

V := E u" + P , (23) 
where e"^ — (dr)"' = d°-i, Ca e° = 1, Co u°- = 0, such that 

VaV = e = - E^ + P'^ , -{Vau'') = E, P={e + E^y/^. (24) 

As spans a radial direction, P > 0. By spatial homogeneity and isotropyQ for a congruence of radial 
normalised geodesies, starting off isotropically from r = 0<^u = (so£' |t,=o = const for all of them). 



= BaE = BaP . (25) 



To determine E, note that 



Sv 



= -VVbiVaU'') = -V^iVbUa)V' 

= -^ehabVV' = -\e>[e+{Vau''f] . (26) 



^See section below. 

®Even though e is not invariantly defined, the 1 + 3 covariant discussion of LRS perfect fluid spacetime geometries given 
in Ref. j?] still applies. As such, is the Fermi-transported (along «") unit tangent of a geodesic and shcarfree spacelike 
congruence. Furthermore, in the given context also its spatial rotation vanishes. 



3 THE GEODESICS 



5 



Thus, we need to solve 



dt 

dv 



V = E = -{VaU") , 



1 dE 



E dE _ I da dt 
(e + £;2) dv ^ a dt dv 



1 d , 



(e + E^) dv 
[e + E^^ 



1 da 
a dt 



-In 

dv 



Integrating, we obtain 

Now solving for E, 
which implies 



(e + Ei) \aJ 



dt 

dv 



y" === E{a) ^ ± 



\ a / 



1/2 



with a "+" for /iiture-directed vectors and a "— " for past-directed ones. Also, with Eq. (|7|), 



^ V habV' = e + E^ = a"^ {V^y 



so 



dr _ _ Pja) 
dv a 



e + E^a) 



1/2 



which, for later reference, can also be cast into the form 

d£:^adr^{e + Eiy/^(^^'^dv 
the definition coming from Eq. (Q). Hence, 



E{a) 



± 



a2 E{a) 



and so 



3.1 Timelike 



dt 
dr 



dt dt I dv 

d^ " d^ ~ ^ [a-2 (e + ^;2(a)) ]i/2 ao{e + E^y/^ 

-1 1/2 



±a(t) 



1 - e 



"0 



■.= ±aoie + Eiy/^ 



(27) 
(28) 
(29) 
(30) 

(31) 

(32) 
(33) 

(34) 

(35) 
(36) 



For timelike vector fields, e = — 1. If we have initially parallel to u'^, then Eq = 1, and so dt/dv — 1 
and dr/dv = 0, confirming that then remains parallel to (which is geodesic). Otherwise, for 
future-directed timelike geodesies that have a non-zero initial hyperbolic angle of tilt with u°' (such 
that Eq > 1), the following relations apply: 



dt 
dv 



1 + iEi 



1) P 
V a 



1/2 



dr 
dv 



{El - 1)V2 ( 



(37) 



3.2 Spacelike 



For spacelike vector fields, e = -I- 1. Setting _Eo = means starting off orthogonally^ but these geodesies do 
not remain orthogonal to the flow lines, and so do not remain within the spacelike 3-surfaces {t — const}. 
Indeed, from Eqs. (p6|) and (P) 



-(i4u")|p = o, e|p>o 



dt 



-(Ku") = E <Q nearby — < 

dv 



(38) 
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showing that the geodesic, nearby spacehke 3-surfaces bend down (into the past) relative to the spaceUke 
3-surfaces {t = const}. In this case ao ^ ±aQ, and 



So, with dr/dv = (ag/a^), we find 



dt 
dv 



dt 
dr 



V a 



-,1/2 



ao 



-,1/2 



(39) 



The geodesic 3-surfaces give the best shcing of a spacetime in order to approximate Newtonian theory in 
a general spacetime — see the discussion by EhlersS — and have been studied in the FLRW context 
by Rindler Page and EUis and Matravers 

The simplest dynamical case is the spatially flat Einstein-de Sitter model, which has (pressure-free) 
incoherent matter as a source, and A = 0. Here, the length scale factor takes the functional form 
a{t) = ao [^Hot]'^/'^, where Hq is the value of the Hubble parameter H := (1/a) (da/dt) at time tg. 
Hence, we obtain ( note that t < to) 



dy 



ao (3i7o/2)2/3 ^ yV3 [ 1 ~ (3Ho/2)4/3 yVS ]l/2 



(40) 



leading to an elliptic integral which gives the value of r at time t, starting off orthogonally at r = and 
time ^0 = I -f^o"^- 



3.3 Null 

In the case of null congruences, e = 0. Then it follows for the past-directed case that 



dt 

dv 



1/2 



= -\Eo\{^ 



and, as dr/dv = \Eo\ (ao/a^), 



dt 

d-r--^^'^ 



(41) 



Alternatively, we can use the fact that ■= ci{t) Ua is a conformal Killing vector field: VaCb = d{t) gab- 
Thus, for any geodesic vector field 



fc''Vfc(eafc") = a{t) k'' gab k' = d{t) ka k'' 
and in the particular case that is null: 

fc^fc" = ^ f„, fc" = a{t)uak"- = const 

^ {kau'^m = 



iak'' = a{t)Uak'' 

const 



Relating this to the redshift, z, defined by 



(1 + ^) 



{kaW") 



a{t) 



E, 



{kbub)o Eo 
for past-directed radial null geodesies it follows that 

T 



a{te) 



k" = 



a{t) 



-1, ^TY' 0' 
a{t) 



- = k' = E=-{l + z) 



(42) 



(43) 



(44) 



(45) 



where we have set i?o = — 1 by choice of the affine parameter v. 
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4 The geodesic deviation equation 

4.1 The deviation vectors 

The basic equations relating the geodesic vector V"" and orthogonal deviation vector 7]°" have been given 
above, see Eqs. (|l|) - (H)- We now restrict the deviation vector further. 

4.1.1 The screen space 

When y is not parallel to w", the vector 77" lies in the screen space of m° (i.e., the spacelike 2-surface 
orthogonal to both, u°- and V"-) iff, additionally to (rjaV^) = 0, r]°- also lies in the rest 3-space of u", i.e., 
iVau"') — 0. We can choose this to be true initially; will it be maintained along the integral curves of any 
geodesic vector field V"'? With Eq. (||), we have 

^-^f^ ^^QKb Tf V"" + Ua rt^bV' = i e /lafe 1/"! + ^f^b^VaU'^) , (46) 

SO 

%!^=^^Vfc(Kzi'^), (47) 
ov 

which will be zero, if 77''Vb(ya'u°) = 0, and this will be true for the congruences we consider (cf. Eq. 
( p5| ) ). Propagation of condition ( ^7| ) along the integral curves of then confirms its preservation. This 
can be seen as follows. The fact that y° and rf" commute, Eq. (||), gives rise to the relation 

= [ u\^,V'') [ybV") + V u'V.VbV" - u^V^rf) (ybV^) - u'^cybV ] , (48) 



which is used to eliminate the respective terms in the "dot" -derivative of condition (47). Hence, with Eq. 

(D, 

[V'Vb{VaU'')-v''^biVaU'')y^lhab[e7^^''V''^ ]-=0, (49) 

which vanishes because hab is symmetric in its indices. So, the consistent solution to these equations is 

(77au'') = 0, (7y,V^°) = 0, D,(nu') = 0; (50) 

i.e., ?7° starts and remains within the rest 3-spaces of u"", and it also remains orthogonal to V^, which 
has a constant scalar product with u° in these rest 3-spaces. From now on we will assume these relations 
hold. 

4.1.2 The force term 

The "force term" (cf., e.g., Pirani jl^ ) for the general GDE (|^) for geodesic congruences of either 
timelike, null or spacelike causal character, specialised to the FLRW case, can now be evaluated from 



Eqs. (|17|) and (|50|) to yield 

R%cdV'rj^V'' ^[e^{fi + A) + ^ {fi + p) E^ ] rj^ (51) 

where, as before, — (Vjju") — E. Note that this force term is proportional to 77" itself, i.e., according to 
the GDE (|^) only the magnitude 77 will change along a geodesic, while its spatial orientation will remain 
fixed|^. Consequently, the GDE (^ reduces to give just a single differential relation for the scalar quantity 
•q. This reflects the spatial isotropy of the Riemann curvature tensor about every point in the present 
situation; anisotropic effects as induced, e.g., by non-zero electric Weyl curvature, Eab, or shear viscosity, 
TTafc, are not involved. 

We deal, now, with three cases: the GDE for a fundamental observer, for past-directed geodesic null 
congruences, and for other families of geodesies. 

^Also, Eq. fclh has no component proportional to u", confirming the consistency of the above screen space analysis. 
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4.2 Geodesic deviation for a fundamental observer 

Case 1: y = u° for the central geodesic. In this case the afhne parameter coincides with the proper 
time of the central fundamental observer, i.e., v = t. From Eq. (pi]), with e = — 1 and E = 1, 



RabcdU^Tl'u'^=\{iJi + ip)ria-\AVa- (52) 

Let the deviation vector be 77" — le'^ , Cac"- — \, Ca = 0, such that it connects neighbouring flow lines 
in the radial direction. Then 6e°'/St — u^Vbe'^ = (as there is no shear or vorticity!), i.e., a basis is used 
which is parallelly propagated along u"^, and Eq. (||) gives 

^^.i^^ + 3p)£+lAi, (53) 

which is the Raychaudhuri equation [ |T6| . However, this equation applies to both comoving matter of 
active gravitational mass density (/i + 3p), and to test matter that is not comoving. On the basis of 
this relation, it is clear that for positive active gravitational mass density and non-negative cosmological 
constant^ all families of past- and future-directed timelike geodesies will experience focusing, provided 
(/i -I- 3p) > 2 A, and so gives rise to the standard singularity theorems ( see, e.g., Refs. ||, ^ ). 

4.2.1 Comoving matter 

For comoving matter, V"' = u"' ^ \ Eq \ — 1 ^ \ E \ ~ 1 ioi the whole family of geodesies. Then, set £ = a 
and multiply by da/dt to get 

da d^a . . ^ . da ^ da ,^ 

Using the conservation equation for comoving matter, 

du 3 da . . dina"^) , „ ^ da 

one finds the familiar Friedmann equation 



^'^\ 1 ,. „2 1 A „2 



2 

I ^a'^ — ^ Aa^ — — k , k — const , (56) 

giving the usual time evolution of a(t) for a given equation of state. In terms of invariants, 

1 da\^ , , . k 

which is just the trace of the Gaufi equation, Eq. (El]), if we identify 



K^^ (58) 

as the constant curvature of the spacelike 3-surfaces {t = const}. Hence, we recover the standard 
dynamical equations for the FLRW models from the GDE. As usual, whenever K is non-zero, by rescaling 
a(t) by a constant the dimensionless quantity k can be normalised to ± 1, which is then the curvature of 
the 3-spaces of maximal symmetry with metric /^j/ dx^ dx'^ ( cf. Eq. (|^) ) . 

If one considers a non-interacting mixture of both incoherent matter and radiation, one has 



3i/o2an„(^) +iHlnr,[^) , p^Hlnro{^) . (59) 



Then, evaluating Eq. (^7|) at i = to shows that 

k 

^0 - i(A*mo +/iro) - |A= -— 4^ Hl{nra„+^ro+^Ko-'^)= Kq , (60) 



'If A < 0, for + 3p) > there will be focusing anyway. 
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where 

Ko-.^^, (61) 

and, as familiar, fli^ denotes dimensionless cosmological density parameters fi; :— /ii/(3i?^) at t = to; 
r^A ■— A/(37?^) defines an analogous quantity for the cosmological constant. Similarly, evaluating the 
Raychaudhuri equation (p3) at t = to gives 



1 /a 



go = 



\a 



= \ ( f^mo + 2 f^ro ^ 2 ilA^ ) ~ i rim^ - JIaj, , (62) 



the to value of the dimensionless cosmological deceleration parameter q :— — {ad a/dt )/{da/dt) . These 



results will be useful in deriving the observational relations for null data (see section 4.3) 



4.2.2 Non-comoving matter 

For isotropically distributed test matter moving with other 4- velocities about the fundamental observers, 
i.e., y = z;" I i?o I > 1, except for the central curve of the congruence w° which coincides with u", we 
need to obtain other solutions to the GDE for timelike curves, evaluated along this central fundamental 
world line (where again proper time t is the same as the preferred affine parameter v, and also here the 
deviation vectors have radial orientation). There are two ways to do this. 

One way is to fully specify the matter source in the equations of the previous discussion on the 
comoving matter case, solve these equations to obtain the source term in the GDE (p3|), and then solve 
the GDE to obtain its general solution (with two arbitrary constants). In the case of the de Sitter 
universe, we have = ii=p, A^d^so Eq. ( |53| ) becomes 

d'^e 

Q^ — -lKl, (63) 



and the solution is 



„, ^ \ Ci cosh(at) + C2 sinh(at) A > , . 

^ ' ~ \ Ci cos{at) + C2 sm{at) A<0 ' ^ ' 

with a := (| | A|)^/^ and Ci and C2 integration constants carrying the dimension of i{t). This shows 
the deviation for arbitrary (i.e., independent of | i?o | > 1) timelike geodesies in the de Sitter (A > 0) and 
anti-de Sitter (A < 0) cases. 

When dynamical matter is present, life is more complex. Defining a dimensionless conformal time 
variable r by dt/dr a =^ d'^t/dr^ — da/dr, for a matter source according to Eq. ( |59| ) the Friedmann 
equation (|^ yields 

da 
d^ 

This can easily be solved when A = 0, for given value of the spatial curvature parameter k. It follows 



[^Aa'^-ka'^ + alH^ n^,a + af, ilr, 1^/^ . (65) 



that the GDE for timclike congruences, Eq. (53), can be rewritten as 



where a = a(r), and da/dr is determined through Eq. (|65|). Unfortunately, this linear homogeneous 
second-order ordinary differential equation is very complicated, except for the de Sitter universe ( where 
= flmg = flro, A^O), which we already considered. 

To provide a simple example with dynamical matter, we fall back onto the Einstein-de Sitter model, 
where A = 0, fc = 0, flrg = fimo = 1. In dimensionless conformal time, the length scale factor is 
a(T) — i Oq Hq t^, and so the solution to Eq. (|66|), which then reduces to 

„ d^i 2 d£ 2 , 



Or, equivalently, (/^ + p) = (/i + 3p) = — 2 fi = const, giving an effective cosmological constant. 
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is given by 



£{t) = Ci r + C2 r2 ; 



(68) 



again, the integration constants Ci and C2 carrying the dimension of £{t) . Fixing initial conditions such 
as to describe a set of test particles isotropically emanating from the central reference geodesic at 77 = 770 1 
one has C2 = — Ci/770. 

Another way to obtain solutions to the timelike GDE ( |53| ) is to use the first integral which relates 
different solutions to the GDE along a central reference geodesic 70 (which is common to both con- 
gruences, and on which the affine parameters coincide and are equal to the preferred time coordinate, 
i.e., V \-yg = t). Let rji relate to the fundamental family of world lines and 772 to another family. Then 
rji = a{t), and as dt/dv = —{vau"-) = E takes the value £^ = 1 on the central reference geodesic, 
{I / rji) (drji / dv) = H ^ ^Q. Considering parallel (radial) deviation vectors for the two families, we 
obtain for their magnitudes 



dr]2 dr]i 



^ ^'^2 const 



In terms of initial data at time t — to, 



(69) 



di]2 
dt 



m I to Ho 



const 
ao 



(70) 



which leads to 



and so 



Then, integration yields 



dr]2 
dt 



1 da{t) uq 
a(t) dt ~ a(ij 



dr]2 
dt 



m I to Ho 



ta 



d 


m 






di 


_a{t) 




a^{t) 


m Ito ( 


'a{t)\ 




dr]2 




+ 


dt 



dr]2 
dt 



m Ito Ho 



.(t) 



a2(y) 



dy 



(71) 



(72) 



(73) 



For the Einstein-de Sitter example, which we referred to before, a{t) — ao {t/to)^^"^ (as Ho = | 
and we find 



m{t) = m Ito ( — 



2/3 



dr]2 
dt 



m kn ^0 ^ 



(74) 



Special cases: 

A: Suppose 7^2 = at i = to (matter flowing out isotropically from the central line at that instant), then 



dm 

dt 



t^/3t2/3(^^-l/3_^-l/3) ^ 



(75) 



to 



giving the radial motion of free particles relative to the fundamental observers, that start off by diverging 
from them. The graph of Eq. (|75| ) was plotted in Fig |l|. 



B: Suppose drj2/dt = at i — tp (matter released from rest at that instant, hence, not comoving with 
the expanding fundamental matter), then 



V2{t) = V2 



1/3 



1/3 



(76) 



gives their radial motion relative to the fundamental observers. The graph of Eq. ( [Tq ) was plotted in 
Fig. I 
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Figure 1: Plot of the deviation vector magnitude r]2{t) according to Eq. ([75|). The parameter values 
chosen are Hq = 60 km/s/Mpc, i.e., to = 0.01 {Mpc/km)s, and dr]2/dt |to = 0.1. 
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Figure 2: Plot of the deviation vector magnitude r]2{t) according to Eq. ([Tq). The parameter values 
chosen are Hq = 60 km/ s/Mpc, i.e., to = 0.01 {Mpc/km)s, and 772 \to — 1 unit length. 

C: Suppose drj2/dt = rj2 \to Hq at t — to (matter initially comoving with the expanding fundamental 
matter), then the matter continues to move as the fundamental observers, i.e., ri2{t) = 772 |to {t/taY'^^ ■ 



Generically, the first integral (|6|), applied to this timelike case, relates the out-going and m-coming 
geodesies that link two (timelike separated) points O and P, on fixing boundary conditions for the first 
integral: namely it relates the positions and velocities of each congruence at O to those at P. Apart from 
the cases just considered, the other one that arises naturally is if particles 1 are at rest at O and coincide 
at P, whereas for particles 2 the situation is the converse: they are at rest at P but coincide at O. Then 



I dri2 



= V2 \P 



dm 

dt 



(77) 



This relates the positions and velocities at O and P, showing that if both distances are the same (in 
absolute, not comoving terms), then the velocities will be the same. 
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4.3 Past directed null vector fields 

Case 2: V = fc", ka fc" = 0, A;'' < 0. Equation 1^ now gives 



(78) 



so writing 77° — rje"", Cae" = 1, = ea w° — Cak"", and using a parallelly propagated and aligned basis, 
Se''/6v = k^Vbe'' = 0, we find from (|), 



(79) 



Again, in line with the timelike case of Eq. ( p3| ) , all families of past-directed (and future-directed) null 
geodesies experience focusing, provided {iJ.+p) > (while the sign of A has no influence). Equation ( [79| ) is 
easily solved in the case of the de Sitter universe, where {^+p) — 0, and the solution is 'ri{v) = Ci V + C2, 
equivalent to the (flat) Minkowski spacetime case. For null rays diverging from the origin, C2 = 0, and 
we have the same angular size-distance relation as in flat space (provided we measure distance in terms 
of the affine parameter v). 

When dynamical matter is present, we need to express the quantities contained in Eq. (|7^) in terms 
of the (non-affine parameter) redshift z, defined in Eq. (ji^). A standard collection of mathematical 
formulae ||l| gives for the derivative operator of Eq. ( [79| ) the expression 



From Eq. (Q) we know that 

(1 + ^) = 

hence, (in the past-directed case). 



a 



E 



dz 



(1 + ^) 



dz'^ dz 



da 
a 



dE 
~E 



1 da 



1 da 



which leads to 



dz= {l + z)- —dv= {1 + z)-—Edv 
a dv a dt 



dv 



EoH{\ + zfdv 



(80) 



(81) 



(82) 



(83) 



dz EoH{l + z)^ ' 

The Hubble parameter is to be determined via the Friedmann equation, Eq. (|57|), from which one obtains 



\^i+\A + H^{l-Qo- f^Ao) (1 + 



By use of the Raychaudhuri equation, Eq. (p3|), one finds, furthermore, 
d'^v 3 . . 1 , „ . 1 



1 



(/i -I- 3p) - 



dz^ EoH{l + z)3 ' 18i?2 ^'^ ^' 9iy2 
So, altogether, the null GDE, Eq. ([tq]), can be expressed in the new form 



A 







d'^T] 



[1 



1 



^ K T drj 1 1 



(84) 



(85) 



(86) 



dz^ (1 + ^) 18i?2 ' 9i?2 " ' dz ' 2(1 + z)2 i/2 

If we consider again the non-interacting mixture of incoherent matter and radiation, we have 

fX^3H^nmail + z)''+3H^nrJl + z)\ p ^ Hf, ^lr,{l + zf . (87) 

Then, from Eq. (jsj), for A = the Hubble parameter evaluates to 

H^^Hl{i + n„,o z + z (2 + z) ) (1 + zf , 
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and Eq. (86) assumes the form 

= 



6 + (1 + 7z) + 2nro{l + 8z + 4z^) d77 
2 {I + z) { \ + n.ni„z + nroz{2 + z) ) dz 



in 



mo 



4a.o(i + ^) 



2 (1 + 2) ( 1 + + Q.roZ{2 + z) ) 



(89) 



When only incoherent matter is present (the dust case), then Vl^o = 0, while a sole incoherent radiation 

0. The popular spatially flat FLRW case is contained for fio = ^^mo + ^^ro = 1- 



matter source has ilr, 



The general solution to this linear homogeneous second-order ordinary differential equation is given 



by 



ri{z) 



(1 



[ Ci ( 2 - a„o - 2r!.„ + n,^oZ) + C2 ( 1 + r!,„„z + nroz{2 + z) f'^ 



(90) 



which we obtained with support from some computer algebra packages. The integration constants Ci 
and C2 carry the dimension of i]{z). With this explicit form for the deviation vector of a (past-directed) 
geodesic null congruence at our hands, we are, now, in a position to easily infer an expression for the 
observer area distance, ^0(2), originally derived by Mattig for the dust case (ri^o = 0), which is of 
considerable astronomical importance (see, e.g., Refs. |l^ and Q ). Using d/dt = Eq^ (1 + z)^^ d/dv 
H {1 + z)d/dz (d. Eqs. ~ ~ 
ri{z = 0) = 0, its definitior 



and (p3|) ) and choosing the integration constants in Eq. (M) such that 



dAo{z) 




viz') 


dQo 




dr,{z')/de\^,^. 



yields 



roiz) := 



ro(z) = [ 2nrno - ( 2 - Qmo ^ ^^ro ) ( ^m„ + 2^!,^ ) 



{1 + zy 



- ( 2 - n„,o - 2nr„ ){! + n^„z + r!,„z(2 + z) fl^ ] , 



(91) 



giving the observer area distance as a function of the redshift z in units of the present-day Hubble radius 
for an arbitrary non- interacting mixture of matter and radiation (and containing as a special case 
the Mattig formula when il,.;, =0). The graph of Eq. ( |9l| ) was plotted in Fig. ^. 

The formula (|9^) is equivalent to the one stated earlier by Matravers and Aziz but — unlike 
the usual calculations — is obtained in a uniform way from the null GDE ( irrespective of the intrinsic 
curvature of the spacelike 3-surfaces \t = const} ). In the usual approach, three separate calculations are 
needed (one for each value of fc) , and it is a matter of some amazement that they all fit the same formula 
in the end. In the present approach, one integration is needed, leading to one formula — a considerable 
increase in clarity. 

The first integral relation can be investigated analogously to the timelike case above. Consider null 
rays diverging from the observer at O and arriving at the source S*, with deviation vector r\\^ and null 
rays diverging from the source S and arriving at the observer O, with deviation vector r/2. The first 
integral is the same as before, but now we need to convert (for past-directed null rays) from the affine 
parameter v to i according to Eq. (^, with ao/a = (1 -I- z). One obtains 



f]i \o 



dm 
di 



m \s 



dm 
di 



(1 + z) , 



(92) 



where the terms drj/di are the angles subtended by the pairs of null rays corresponding to the deviation 
vectors. Expressed in terms of angular diameter distances, ro and rs, defined by 



m \s ■= ro 



dm 

d£ 



m \o 



rs 



o 



dm 

d£ 



(93) 



-'dQo here denotes an infinitesimal solid angle rather than a change in density parameter. 
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Figure 3: Plot of the observer area distance ro{z) according to Eq. (91), in units of T/q ^. The parameter 
values chosen are Hq — 80 km/ s/Mpc, ilmo — 0.2 and = 0.1. 



(which, for FLRW geometry, are the same as area distances), we find the familiar null reciprocity theorem 
for FLRW models |: 

rs = ro (1 + z) . (94) 

This underlies the equivalence (up to redshift factors) of area distance and luminosity distance, and the 
fact that measured radiation intensity is independent of area distance, depending only on redshift ( see 
Ref. Q for a more detailed discussion). These features are fundamental in analysing observations of 
distant sources and measurements of the cosmic microwave background radiation. 



4.4 Generic geodesic vector fields 

Case 3: Va V"" = e, not parallel to m", nor null. The force term in the generic case is provided by Eq. 
(^). Writing 77° — £e°-, ea e° = 1, = CaV"' — Cau"--, and employing a parallelly propagated and aHgned 
basis, Se'^/Sv = V^Vbe" = 0, we find from Eq. (|), 

= -6i(M + A)^-i(M+p)£;2£, (95) 

giving the spatial orthogonal separation of these geodesies within the 2-D screen space as they spread 
out in spacetime. 



4.4.1 Orthogonal spacelike geodesies 

A particular case is the spatial geodesies that start off orthogonal to u° (so Eq = 0, which implies that 
the corresponding geodesies are indeed spacelike), but then bend down towards the past thereafter (see 
the discussion in section |3^ ) . The above equation applies with e = 1. The simplest case is a de Sitter 
universe where = /i p, and then the solution for all \Eq\>Q, i.e., all spacelike geodesies is 

Ptn^ = / ^) + sin(a v) A > . . 

^ ' \ Ci cosh(aw) + C2 sinh(at;) A<0 ' ^ > 

with a := (^ I A 1)^/^ (note this is just the exact converse to the timelike case of Eq. ( |6^ ) above). 

In the case of non-zero dynamical matter, however, /i and p are not constants along the initially 
orthogonal geodesies, as these geodesies do not remain within a spacelike 3-surface {t — const}; we have 
to find ^[t{v)] or /i[t(r)] from the geodesic equation. However, near the starting point t; = at we 
have (for A = 0) 

i^o^, (97) 



to 
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giving the solution near this origin, on carrying out a first-order expansion, by 

d£ 



i(v) = £q cosfwo v) . 

dv 



sin(t^o«), wo (i/^o)'/' • (98) 



This is always convergent for normal matter, irrespective of the intrinsic curvature of the particular 
spacelike 3-surface {to — const} considered. However, as soon as the distance is appreciable, the geodesies 
will have bent down and lie below the initial 3-surface {to = const}, where the density of matter will be 
higher and the curvature greater. Thus, the geodesies will tend to converge even more strongly. 

4.4.2 Geodesies in the orthogonal spacelike 3-surfaces 

This is to be contrasted with geodesic congruences within the spacelike 3-surfaces {t = const} orthogonal 
to u", which are 3-spaces of maximal symmetry. In contrast to Eq. (|]), these geodesies satisfy the 3-D 
equations 

V — , K V'" = 1 , K = 0, = V^V>bV'' . (99) 
dv 

From Eq. ( |2^ ) , the force term for the resulting 3-D spatial GDEp] takes the form 

"'RabcdV'l^^V ^\{^i-\Q^ +K)lla^Kl^- , (100) 

where K{t) is the curvature of these 3-spaces (cf. Eq. ( pi] ) ). Consequently, whether geodesies in these 
spacelike sections converge or diverge depends on the sign of K. Setting 77" = 776° where Ca e** = 1 and 
Ca u° = 0, as before we choose a congruence of vectors such that Se"" /5v = y''Df,e° = and the 3-D 
spatial GDEM becomes 

J-'^"^ (101) 

K = K{t) is indeed constant along these spatial geodesies (because they lie within the 3-surfaces {t = 
const}). If i^T > 0, one deals again with the familiar oscillator equation, i.e., two neighbouring spatial 
geodesies will harmonically converge to and diverge from each other as v increases, li K < Q, they will 
exponentially diverge, and if = 0, they diverge linearly. 



Focusing on radial spatial geodesies, the local FLRW coordinates of the spacelike 3-surfaces {t — 
const} arise as follows. We consider a 3-space with metric f^u dx^ dx'^ , and constant dimensionless scalar 
curvature, if non-zero, normalised to A: = ± 1, (cf. Eqs. (^ andjs^) ). Note that the full 3-space metric 
hfj,u{t) at arbitrary time t is just given by hf^^{t) = a? (t) /^t^ .p| Choosing an affine parameter v — r, 
V"" ~ ((^r )" = <^°i is the geodesic unit normal to the 2-surfaces {r — const}, which are 2-spheres of area 
47r/^(r). Thus, it is tangent to the orthogonal coordinate curves x^ = const, x^ = const. A basis of 
deviation vectors in the 2-D screen space is given by 771° ~ (5 "2 and 772° = ^% (these commute with the 
geodesic vector V"' = S^, because each of these is a coordinate basis vector). Employing an orthonormal 
basis with components (ei)" = 5\, (£2)" = f~^{f) S°'2, (ea)" = f~'^{r) (sin 6')"^ 6%, parallelly propagated 



along the radial geodesies Eq. ( |lOlD yields 



the second relation following because relative to the orthonormal basis, 771" = /(r) ^"2 and 7/2° = 
/(r) sin 9 6% (apply the first equation to either vector to get the second). Then the solution we want 
corresponds to that solution for which 77(r = 0) = 0; we find 

sin r fc = + 1 

fir) ={ r k = , (103) 

sinh r fc = — 1 



iiDetermined by Eqs. (g), (U), andW. 

'^^That is, the 3-D version of Eq. (H) that applies in these 3-spaces. 



'^^When a(t) is of unit magnitude, say at time t = i, then ff^v is equal to the metric h^v(t) on the 3-surface {i = const}, 
except for a dimensional unit factor, and similarly for k and Kit). 
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showing how the GDE withm the spacehke 3-surfaces {t — const} determines the function f{r) in Eq. 
(0). The corresponding solutions with drj/dr = at r = exhibit precisely how Euclid's parallel postulate 
breaks down for these curved 3-space sections, according to the spatial curvature. 

In this context it is of interest to remark that the Lorentz-invariant|^ de Sitter spacetime geometry, 
which is the case = ^ — p, A > 0, can be sliced by spacelike 3-surfaces {t = const} of either constant 
positive, zero, or negative intrinsic curvature (cf. Ref. ), depending on the sign of the sum 3K = 
— ^ 6^ + A (see Eq. ( ^ ) ). For anti-de Sitter (A < 0) only the negative curvature case applies. The 
different FLRW forms of the de Sitter spacetime metric follow from arguments essentially identical to that 
just given for the 3-space metric, because it is a 4-space of constant curvature, i.e., maximal symmetry 
( and the argument applies also to the 2-sphere, leading to the form of the terms in the last bracket in Eq. 

). In each case, the GDE, together with the constant curvature condition (20), leads to the harmonic 
equation (102). 



Similarly to the null case, the 3-D geometrical reciprocity theorem can be stated as 

dri2 



m \o 



dr 



m \p 



dm 

dr 



(104) 



showing how geodesies diverging about a central geodesic from P to O at an angle ao reach a separation 
d at O, and geodesies diverging from O at the same angle will reach the same distance apart at P 
(irrespective of the spatial curvature which is constant). Corresponding statements hold for the families 
of geodesies that diverge from P and O, and end up parallel at O and P, respectively. 



5 Conclusion 

One way of solving the EFE is to treat them as algebraic equations relating Rabcd to Rat and Cabcd, then 
solving the GDE (which characterises relative acceleration due to spacetime curvature) to determine both 
the spacetime geometry and its properties. In the case of a FLRW model, this can be carried out explicitly, 
as shown above: integrating the GDE (cf. Eqs. (|5^), ( ^6|) and (|95|) ) allows complete characterisation 
of all interesting geometrical features of the exact FLRW geometry in an elegant manner — determining 
the timelike evolution, spacelike geometry, and null ray properties, which in turn determine the basic 



observational properties. The Newtonian analogue of some of this has been given by Tipler |27|, |28 . 

An interesting project is to extend this calculation to perturbed FLRW models in order to work out 
the effects of linear anisotropics on the present results as regards all three causal cases (timelike, spacelike, 
null). This would allow investigation of both dynamical and observational features of such models, for 
example examining aspects of gravitational lensing theory pO| . 
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